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This paper considers the conductance through a dirty junction between a normal conductor 
and a superconductor with an insulating barrier. It is shown that for large barrier resistances there 
is a relative enhancement of the conductance near zero voltage, whereas for low barrier resistance 
this anomaly appears at finite voltage. 
PACS numbers: 74.80.Fp, 74.50. +r 

Recent experiments on contacts between semiconductors and superconductors have generated a lot of interest 
] (see e.g. and references therein). Besides possible technical application of these devices, one is interested in 

understanding how electrical conductance below the superconducting gap arises. While there is little doubt that 
the basic process is Andreev reflection [|j , the role of impurities is only being appreciated gradually recently. The 
^ , interplay between impurity scattering and the potential scattering at the Schottky barrier between the semiconductor 
O and superconductor is believed to be the physics behind the zero bias anomaly (ZBA) observed in many 

. experiments (e.g. [Q-§)) 

To understand the transport in these normal-insulator-superconductor structures two major theoretical tools 
have been employed. In Ref ||l3,0 the conductance is expressed in terms of the transmission matrices across the 
structure. Unfortunately quite often the necessary averages are not available from random matrix theory and thus 
numerical calculations to obtain the transmission matrices have to be done for a large number of "samples" and then 
' averages over the samples taken. An alternative approach |]7|-|l0|| is to use the quasiclassical green's function where 
the impurity average has already been made. However, in this approach the impurity scattering is only included in 
■ the self-consistent Born approximation. Moreover these references confined themselves to the cases of high potential 
. barriers (and often low energies) so that analytical progress can be made, although the basic equations are applicable 
^-H ' more generally so long as we are in the dirty limit. In this paper we shall solve these equations numerically in order 
to obtain a better understanding of the physics as well as a detail comparison with the results from the scattering 
matrix approaches. | |r^ , [T^ 

We shall then consider a (quasi-one-dimensional) normal conductor (" wire" ) situated between a normal [x — 
0) and a superconducting {x — L) reservoirs. The reservoirs are assumed to be at equilibrium at potential V and 
I ■ respectively. A potential barrier is located at where < < L. We shall confine ourselves to the dirty limit 
' where all the energy scales are smaller than the impurity scattering rate. We also require that the mean free path I 
,: to be much less than L. (Note that experiments cited above may not satisfy this condition). We shall assume that 
there is no pairing interaction in our conducting "wire" and we shall consider only ordinary and magnetic impurity 
scatterings. To study the transport through the system we shall solve the Usadel equation 



On 
On 



(N 



0\ 



o 



[era - dspm^g] + —df,{gd^g) (1) 

TT 

together with the normalization condition g^ = — tt^I governing the angular averaged Green's function 

g^ 

and g^'^'^ are the retarded, advanced, and Keldysh martix Green's functions in particle-hole and spin space. Here 
e is the energy and ^fp^'^ — ^T^g^'^'^T^ is the self-energy for magnetic scattering with the scattering rate 7. The 
equation connecting g on the two sides of an interface has been derived by Kuprianov and Lukichev | [l5| ] , 

aSgd^g = ^[g{xb-),gixb+)] (2) 

a is the conductivity of the wire, S the area, and is the conductance through of the barrier in the normal state. 
For our purpose g^ can be chosen to have only T3 and <T2Ti components: g^ = —477(93X3 — /i(72Ti). g^ is related to 
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by = T3(^^) V3. Since g| + = 1 from the normalization condition, we can write 53 = cos6' and /i = sinff 
where 6 = 6r + 0i is in general complex. From eq(|^) 9 satisfies 

2i{e + ncose)sme + Ddle ^ (3) 

The boundary condition eq(||) reduces to 

{aS)d^9 = G^in([0]f,) (4) 

here [9]i, is the discontinuity of 6 across the interface, i.e., [9]i, = 0{xb+) — 0(.Tf,_). 6 at x = and x — L are assumed 
to attain their equilibrium values of the reservoirs. Thus at x = Q 9.^ = 9i = Q and ai x = L, 9^ = ^,9i = jln^^f, 
i.e., tanh6l, = ^ for e < A. ( and 9r = 0, 6I, = |lnf^, for e > A ). 

We can define a distribution function h by g^ — g^h — hg^ where h can be chosen diagonal h — /iqTo + /13T3. 
For our purposes it is sufficient to solve for hs, which satisfies (from eq(|l|)) 

d,icosh'9Ah3) = (5) 

in the bulk together with the boundary condition (from eq(^)) 

{aS)cosh'^9idxh3 = cos([6'r]fc)cosh6'i(a;f,+ )cosh6'i(a::b_)[ft.3]b (6) 

Here [hsjb is the discontinuity of /13 across the barrier. For a potential V applied to the normal reservoir at a; = 
and zero potential at the superconducting reservoir, we have then /13(e) = _ q q^^j^ h^ie) — 

at X = i. Here n{e) = tanh2^. We shall confine ourselves to zero temperture, where n(e) becomes a step function. 
Thus /13(e) decreases from 1 at a; = to at x = i if |e| < eV, and is identically zero otherwise. Since the current 
/(V) = Sf J^^^tr[T3{g'^d^g^ + g^d^g"^)] can be written as I{V) = decosh^9^d^hs„ the conductance 

G{y) = ^ at voltage V is given by 

G{V) = -dScosh'9Ah3\,=,v (7) 

It is obvious from the discussion above and eqs(||,^) that /13 behaves like a voltage along a wire with ( position 
dependent) conductivity aS'cosh^^i in series with a resistor with conductance 

Gt{V) = GlcoB{[9r\b)cosM^{xb+)cosM,{xb-)\,=eV (8) 

and can be easily found. /13 of course should not be taken literally as the voltage itself, but the concept of these 
effective conductances will help us understand our results. The conductance through the system can therefore be 
written as 

G{V) = {Ri{V) + Rt{V) + R2{V)r^ (9) 

where 



R^{V) = ^ / ' dxsech^9,{. 
Jo 



x)\e=eV (10) 

and similarly for R2 except the integration limits are from Xb to L and Rt = {Gt{V))^^ . The above equations have 
essentially been derived in ref. 

Before describing our results it is helpful to discuss the physical meanings of 9r and 9i. The usual density 
of states is given by Regs = cos9rCOsh9.i. Notice however this is not the quantity that is directly relevant to the 
effective conductances through the wire or the barrier. As seen from eqs(^0|) and (||) the relevant factors are cosh^^j = 
(Regs)'^ + (Re/i)2 and cos{[9r]b)cosh9i{xb+)cosh9i{xb-) = {Reg3{xb-))(Reg3{xb+)) + {Refi{xb-)){Refi{xb+)) respec- 
tively. Notice also under a unitary transformation between the particle and hole components of opposite spin (i.e., a 
generalized Bogoliubov transformation), g^ transforms into g'^ — e'^~"'^'^^g^e~^^'^^'^^ where x is a real number. If one 
parametrizes g'^ by 9' as before then 9'^ = 9r + X ^-nd 9'^ = 9i. It is thus useful to think of 6*^ as the angle of a vector in 
an imaginary space which I shall refer to as spectral space (see also p^), and cosh9i as its magnitude. The equation 
of motion for 9r and 9i is given by eq(^ and Regs and Re/i are then the two {g- and /-) components of the 
spectral vector or the density of states vector (DOSV). The effective conductance of the conductor at a given point 
and voltage is thus enhanced by the square of the DOSV (at that point and voltage), while the conductance through 
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a barrier is modified by the factor corresponding to the scalar product of the spectral vector on the two sides. Both 
these factors are invariants under the particle-hole transformation. Conductance through a barrier is thus possible 
only when the spectral vectors on the two sides have finite projection onto each other. For a direct contact between 
a normal metal {6r = 0) and a superconductor {9i — for e < A) there is thus no subgap conductance. In the NIS 
structure as in here the finite resistivity of the wire allows 0^ 7^ at Xb- and hence a finite barrier conductance. 

We shall concentrate on the subgap conductance eV^ << A. For a better display of the 9 parameters in Fig. 
^ below, I have chosen to present the results for Xb — 0.92L. The results for the total resistances presented below are 
virtually independent of Xb so long as Xb ~ L. From the above equations it is convenient to define a dimensionless 
parameter gb as the ratio of the normal state conductances through the barrier to that of the wire, i.e. gb = /G^^^^, 
where G^i^e — o'S/L. We shall also express all energies in terms of the natural scale Eo = D/LF'. The gap A is 
chosen to be IQQEd. 

Fig. |l| shows the conductance of our NIS system for intermediate values oi gb- For (^b^l.O we obtain the an 
anomaly at zero bias (ZBA). At higher g^'s an anomaly at finite voltage (FBA) is apparent. The position of this FBA 
is roughly ~ gbEn for these intermediate values, (see however below). This FBA is actually built on top of a ZBA. 
Fig. H shows the effect of magnetic scattering on the FBA for gb — 2.5. We see that the FBA is sensitive to 7, and is 
removed for 7 ~ gbEn (also see, however, below). The ZBA remains after the FBA is destroyed. 

To understand these features, it is helpful to consider the 'resistances' defined in eqs (||)and (p^), as shown in 
Fig. H again for 5^ = 2.5. (l^ The FBA results from the dip at finite voltage for Ri (more precisely Ri + R2, the total 
effective resistance for the wire), whereas the ZBA is from the minimum at zero voltage for Rt- The relatively small 
7 ( ~ gbEo for these intermediate gb) first remove the feature in The minimum at ^ = for Rt is somewhat 
more robust. 

The features for i?i and Rt can in turn be related to the behavior of 9r and 9i. (recall eqs(||) and (|lO|). It is 
also helpful to observe that a natural length scale i.e., the coherence length at energy e, arises in eq (^) ) Fig. 
^ shows these parameters again for gb — 2.5 and for a specific energy e = 2.5. First consider 7 = 0. We see that 9i 
is finite within the wire, but small at the two ends {9i w at a; = L because e << A). This 9i is generated from 
the gradients, as is evident from taking the imaginary part of eq(|^) : at finite e and 0,. there is a contribution to a 
negative curvature of 9i in the form of jysm9rCOsh9i. Fig. || shows the values of 9r and 9i at the immediate left of the 
barrier (xb-) for two values of gf, at 7 = 0. These values are good indicators of the values of 9r and 9i within the wire 
if we keep in mind their typical spatial dependence as in Fig. ^. As the energy increases 9i{xb-) first increases, then 
saturates and shows a slight decrease afterwards. This gives rise to the shape of At larger energies 9r decreases 
faster from the boundary towards a; = 0, thus acquiring a larger slope at the boundary and hence a larger [9r]b (as well 
as a smaller value within the wire which slow the increase of 9i as just mentioned), and gives rise to the increase of 
Rt as a function of voltage. For smaller gbS the resistance is dominated by Rt, where the angle between the spectral 
vectors on the two sides of the barrier are also larger. Thus for smaller gb as a function of voltage Rt increases rapidly 
with voltage and dominates the resistance, and only the ZBA but not the FBA remains. Thus in short the ZBA is 
due to the tilting of the spectral vector at Xb- at low voltages alway from 9j. — due to the proximity effect , and 
the FBA is due to the fact that at finite energies a spatial gradient increases the magnitude of the DOSV within the 
wire. The presence of a small 7 suppresses the negative curvature of 9i and thus 9i itself, thus destroying the FBA. 
7 also increases the curvature of 9r and thus [9r]b, which contributes to the suppression of the ZBA. (see Fig. ^ ) 

An example for the conductance at small gb (0.1) is shown in Fig. ||. jl^ It is interesting to note that the 
curves for different 7's actually cross each other at cV^Ed, i.e., for some "large" voltage that G actually increases 
with increasing pair-breaking. This anomalous dependence has also been found by Volkov and Klapwijk |^ for the 
double barrier structures for large barrier resistances. It is not clear how to understand this physically. 

For very large gb the resistance of the system comes almost entirely from the wire, and the FBA is most 
prominent. The position of the FBA as well as the 7 needed to destroy the FBA seems to saturate for large 7 (at 
5Et> and ^ 'SEd respectively): see Fig. ^ (^f, = 24). Since increasing e tends to increase 9i, for large 7 the 
conductance actually has a minimum at eV = and increases with increasing voltage. 

It is of interest to compare Figs |^J^,^,^ with the results of Marmorkos et al iQ , who have studied the ZBA 
and FBA by employing the scattering matrix method. The qualitative behaviors of the conductance are similar. Our 
FBA seems to be located at a similar value of eV and has similar width |l^. However, important differences have to 
be noted. It is useful first to recall that the normal state conductance of the system is GJ5,j^g/(l -I- g^^) ( e.g, O.SG^^^^, 
for gb — 4.0 and 0.96G'^^j.^ for gb = 24). Hence for gb^4: the conductances at eV = are almost the normal state 
values, and for the voltages shown for these f;f,'s the conductances are always larger than the corresponding normal 
state values except for very small voltages. The conductances at eV — found in ||l^ are smaller than the normal 
state even for perfectly transmitting interface, and only overshoot the normal state value near the FBA for relatively 
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transparent barriers pOt . 

Marmorkos et al and Takane and Ebisawa (TE) had also examined the effect of magnetic field B 
on the conductance at zero voltage. The magnetic field makes the problem intrinsically (at least) two dimensional in 
general and we hope to return to that in the future. However, since the magnetic field is also a pair-breaker it is of 
interest to compare their results with Fig. ^, ^ and|7[ It is obvious from our plots that G{V — 0,7) always decreases 
with increasing 7. The relative change in G (i.e. the change measured in G{V = 0,j)/G{V = 0, 7 = 0) ) decreases 
as the barrier resistance decreases, and vanishes as — > 00 (see Fig. |, and also . In TE the behavior of 
G{V = 0, B) as a function of increasing magnetic field B is similar to the one for G{V ~ 0, 7) just described. However 
in Marmorkos et al although G{V = 0, B) still decreases with B for high potential barriers, it increases with field for 
the transparent interface. It is unclear why the results of Marmorkos et al and TE are different. 

For normal metals if ~ lO^cm/sec and L ~ ^m, then Eu ^ 0.02meV if / ~ lOOOA. To observe the FBA 
the condition > GJ^j^e roughly requires the transmission probability (appropriately angular averaged) through 
the barrier to be j^O.l, which should be easily achievable for a suitable choice of the metal and the superconductor. 
To observe the FBA in the semiconductor-superconductor system seems to be not so easy. Even for a choice of 
semiconductor like InGaAs where the tunneling through the Schottky barrier is enhanced by the small effective mass, 
since the Fermi wavelength on the two sides differ significantly most of the excitations are still reflected. In order to 
get > Gj^j^g smaller / and/or larger L have to be chosen, making Ed much smaller. One should also keep in mind 
that most of the recent experiments on these structures involve semiconductors in the clean limit, where the present 
calculation does not apply. 

This research is supported by the National Science Foundation through the Northwestern University Materials 
Science Center, grant number DMR 91-20521. 
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~ 3.14Sd here). 
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(1979) (See also Y. Takane and H. Ebisawa, J. Phys. Soc. Jpn 61, 2858 (1992); C. W. J. Beenakker Phys. Rev. B 46, 
12841 (1992), and ref ) that the resistance of the wire connecting N and S for \/ is the same as that in the normal 
state in the limit L » L Thus this difference may be due to a finite l/L correction . 
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FIG. 1. The (normalized) conductance versus voltage for gb's listed in the legend. 7 = 




FIG. 2. The (normalized) conductance versus voltage at gi = 2.5 for the values of 7 listed in the legend. 
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FIG. 3. The effective "resistances" Ri and Rt for gt, = 2.5 for some values of 7 listed in the legend. RZire = (GJ^ire) ^• 
i?2 is very small, roughly constant (~ 0.08i?I];j^e) and is not shown here. 
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FIG. 4. The parameters 6r and 9i for gi, = 2.5 and e = 2.5 for some values of 7 listed in the legend. 
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FIG. 5. The parameters Or and 6i at Xb- for gt = 2.5 and gb = 0.5 as a function of energy. 7 = 0. 
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FIG. 6. The (normalized) conductance versus voltage at gb = 0.1 for the values of 7 listed in the legend. 
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FIG. 7. The (normalized) conductance versus voltage at gi, = 24 for the values of 7 listed in the legend. The horizontal 
dotted line is the normal state conductance of the system. 
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